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NONSYMMETRIC GENERIC MATRIX EQUATIONS 



GERALD BOURGEOIS 



Abstract. The n X n matrices A, (B;);,C are indeterminates. Let K = 
Q(E), where E denotes the entries of the previous indeterminates, and let 
K be the algebraic closure of K. We show that the Riccati equation XAX + 
B±X + XB2 + C = n has ( n ) solutions in the unknown X £ A4n(-K") and is 
solvable by radicals over K if and only n = 2. Solving the previous equation is 
equivalent to find the roots of a polynomial whose Galois group over K is S2n ■ 



We show that a polynomial equation g(B\ 



,B k ,X) + B, 



fc+i 



n , where g 



has no constant term, admits a finite number of simple solutions. We give the 
Hilbert dimension of the variety of solutions of the equation X 2 + B\ X = n . 
We consider the unilateral equation X k + Bj,_ 1 X k ~ 1 + • ■ • + B\X + Bo = n 
when the n X n matrices (B;); are commuting indeterminates ; we show that 
every solution X £ M n (K) commutes with the (Bi)i. When n = 2, we prove 
that the equation B±XB 2 X + XB 3 X + X 2 B 4 + B 5 X + B a = 2 admits 16 
isolated solutions in M2(K). Again when n = 2, we present the different 
forms of the set of solutions of a numerical Riccati equation. 



1. Introduction 

S. Gelfand wrote in 2004 (cf. [5]) : "The problem of solving quadratic equations 
for matrices seems to be absolutely classical. It would be natural that such a 
problem should have been at least formulated, or even solved, in the 19th century 
at the latest. Still, I asked many people about this problem, and they directed me 
to various sources, but nowhere could I find even a mention of this problem" . 
In this paper, we deal with generic matrix equations, that is, polynomial equations 
whose coefficients are indeterminates. Questions about generic matrices are solved 
in [I] and [2] or about formal matrices in [T3]. More generally, C. Procesi described, 
in |15j . properties of the algebra of generic matrices. 

Let Q be the field of rational numbers. If M is a n X n matrix, then xm denotes 
its characteristic polynomial, a(M) its spectrum and tr(M) its trace. 

Definition. The n x n matrices {A r = [a r; ij])i<r<fc ore non-commuting indeter- 
minates and the entries of which are commuting indeterminates. We consider the 
quotient field K = Q.((ai- t i t j)ij , • • • , {o>k;i,j)i,j) o-nd its algebraic closure K . Let f be 
a non-zero polynomial over K in fc + 1 non-commuting indeterminates. We consider 
a matrix equation in the form 



(1) 



f(A 1 



, Ak,X) = n in the unknown X = [xi.j] € M n (K). 



i) Assume that the previous equation has a finite positive number of solutions. If 
the entries of each solution can be calculated by radicals over K , then we say that 
Eq (l\) is solvable, else we say that Eq {!]] is non-solvable, 
ii) (cf. [12j ) A solution X$ of Eq (QP is called (geometrically) isolated if there is a 
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2 GERALD BOURGEOIS 

neighborhood of Xq that contains no other solution of the equation. 

df 
A solution Xq of Eq fip is called singular if det(— — (A\, ■ ■ ■ ,Ak, Xq)) = 0. 

oX 

First, we study the nonsymmetric algebraic Riccati equation in X € M. n {K) 

(2) XAX + B X X + XB 2 + C = 0„. 

where the n x n matrices A,B\,B2,C are indeterminates. For the general study of 
Eq ©, in Section 2, there are two parts 
i) One reduces Eq J2J to the equation 

(3) X 2 + BX + C = n . 

where the n x n matrices B,C are indeterminates. In |18j . J.J. Sylvester studied 

Eq©. 

ii) In Theorem [TJ we prove that Eq © has ( ) solutions and that solving Eq ^ 

is equivalent to find the roots of a polynomial of degree ( 2 ^ 1 ) that has S^n as Galois 

group ; thus the Riccati equation is solvable if and only if n — 2. 

In Section 3, we give two instances of specializations of Eq ([2]) such that they 

are solvable until n = 4. In particular, we study in details the generic equation 

X 2 + BX = n . 

In Section 4, we show (Theorem^ that an equation of the form g(A% , • • • , Ak, X) + 

B = 0„, where the polynomial g has no constant term, and the (Ai)i and B are 

indeterminates, has a finite number of solutions in A4 n (K) and these solutions have 

multiplicity 1. 

To obtain explicit solutions, we use in Section 5 and 7, for low dimensions, Maple 

and Magma softwares ; this method permits to calculate the minimal polynomials 

of the entries of the solutions of Eq © and also to study non-unilateral quadratic 

equations. We consider the system of algebraic equations, in the unknowns {xij)ij 

and the parameters (bi,j)i,j, ( c i,j)i,j> associated to Eq (|3|). If n = 2, then we can 

completely solve this system. Yet, when n = 3, 4, we must use specializations of 

B,C. 

Again using formal calculus when n < 3, we have a look at the non- unilateral 

equation 

(4) X 2 + B 1 XB 2 + C = n . 

The software is unable to solve the generic equation ; yet we can test specializations 
of Eq ©. We show that, when n > 2, this equation admits a finite number of 
solutions and is non-solvable. 

In Section 6, we study the unilateral matrix equation X k + Bk-\X k ~ l + ■ ■ - + B\X + 
Bq = n when the n x n matrices (-Bj)i are commuting indeterminates ; we show, 
in Theorem^ that the solutions X € A4„(K) commute with the (Bj)j. This result 
works for commuting Riccati Eq © or commuting Eq (£[]) and does not work for 
the equation 

X 2 + BXB + C = 2 

where the 2x2 matrices B, C are commuting indeterminates. 

In Section 7, Theorem 2] gives, for n = 2, an instance of a generic quadratic matrix 

equation that admits 16 isolated solutions in K , that is, according to the Bezout's 

Theorem (cf. [6]), the maximal number of isolated solutions of a quadratic equation 

in 2 x 2 matrices. 

In Section 8, again when n = 2, we present the different forms of the set of solutions 

of a non-generic Riccati equation. 



nonsymmetric generic matrix equations 3 

2. The nonsymmetric algebraic Riccati equation 

In this section we will see that the complete Eq © is not solvable for n > 3. Note 
that the form of Eq (|2|) is invariant by translation of the unknown. Let X = Y + U 
where U is a constant matrix to be chosen. We obtain Y AY + (UA + B\)Y + 
Y(AU + B 2 ) + UAU + B X XJ + UB 2 + C = 2 . Since A is generic, it is invertible 
and we put U = —A~ 1 B 2 . Thus the study of Eq © is reduced to the study of the 
following : let A,B,C be n x n matrices that are indeterminates. We consider the 
equation in the unknown A G M. n (K) 

(5) XAX + BX + C = 0„. 
Since A is invertible, Eq ([5]) is equivalent to 

(AX) 2 + [ABA- l )(AX) + AC = 0„. 

Finally the study of Eq ([2]) is reduced to the study of Eq ([3]) in the unknown 
X = [x%j] G M. n (K), where B, C are n x n matrices that are indeterminates. 
Yet, for the study of the complexity of solving Riccati equation, it is easier to 
consider Eq (J2J). 

Lemma 1. Let H — \jtij] be a p x p matrix where the (^i,j)i,j are commuting 
indeterminates and K — Q((TTij)i,j). Then \n is irreducible over K and its Galois 
group is S p . 

Proof. Let P be a polynomial of degree p with coefficients in Q that has S p as 
Galois group. We specialize II into IIo so that IIo is the companion matrix of P. 
According to [20j ch. 8.10], S p is a subgroup of the Galois group of Yn and we are 
done. □ 

Theorem 1. Let n > 2. Let A, B\, B 2 , C be nxn matrices that are indeterminates. 

Then 

i) Eq 0), in the unknown X, has exactly ( ^) solutions that are in M. n (K((\i)i)) 

where the (\i)i< 2n are the roots of a polynomial that has S 2n as Galois group over 

K. 

ii) Eq (0) is solvable if and only if n = 2. 

Proof, i) Following the idea of J. Potter and B. Anderson (1966), we associate to 
Eq (12) the generic 2n x 2n matrix 

M=(- B2 ~ A 
\ C B x 

If A is a solution of Eq @, then the graph of A is M- invariant (cf. [H]). This 
essential property works because the LHS of Eq @ is a linear function of the inde- 
terminates. More precisely, according to [5], there is a one-to-one correspondence 
between the set of solutions of Eq (J2J) and the set of n-dimcnsional M-invariant 
subspaces E that satisfy the condition 

(6) E is complementary to {0} x K . 

According to Lemma [U xm is irreducible over K and its Galois group is S 2n . Then 
M has 2n distinct eigenvalues (Aj)j that are conjugate over K . Thus we can choose 
an associated basis of eigenvectors B = {/i, • • • , f 2n } G L 2n where L = K((\i)i) is 
a finite extension field of K and such that the (fi)i are coordinatewise conjugate 
polynomials in the (Ai) with entries in K. Any n-dimensional M-invariant subspace 
is spanned by n vectors /j 1 , • ■ • , /j n , where i\ < • •• < i n , chosen in B (cf. [H ch. 

2.1]). Put (f h ,-■■ ,fi„) = [y J ) wner e J = {h, -•• , ««} and Uj, Vj are n x n 

matrices. Then the subspace [/i 17 • • • , /»„] satisfies Condition ^ if and only if Uj 
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is invertible. Assume that the polynomial, in the (Aj), det (t/j) is zero. We consider 
a symmetrization of det (Uj), H cr£ s 2n det(C/ cr ( il ) i ... )CT (j„)) = 0, a polynomial relation 
linking indctcrminatcs in K. We obtain a contradiction if det (t/j) is not identically 
zero, that is, if there is a specialization of M such that the associated det(t/j) is 
not zero. 

We specialize M into Mo, the anticirculant matrix associated to the sequence 
1,2,- •• , In. We can show that Mq has distinct eigenvalues and, for any choice 
of n eigenvectors of Mo, the associated matrix Uj is invertible. Note that every 
matrix in a neighborhood of Mq in M.2n{Q) satisfies the previous properties ; in- 
deed, since Mo has distinct eigenvalues, cr(M) and B are continuous functions of 
M in a neighborhood of Mo. 

Therefore, any choice of n-dimensional M-invariant subspace is convenient and Eq 
© has exactly ( 2 ^) solutions. Moreover, the solution associated to [/j i; • • • , /»„] is 
VjUj , that does not depend on the choice of the ordering of the eigenvectors ; 
we can even replace fi 1 , • • • , fc n with any basis of [fi 1 , • • • ,fi n ]. Thus, if we know 
explicitly the (A-j),, then the calculation of the solutions of Eq @ is reduced to 
algebraic operations in L. 

ii) Since S4 is solvable, the result is clear when n — 2. 

Let n > 3. Since S^n is non-solvable, it remains to show that if the solutions 
(-Xj) .wan} of Eq © are in M„(L) where L D K, then L D K((Xi)i<2n)- Ac- 
cording to [5], a n-dimensional M-invariant subspace, associated to a solution X, is 

5(x) = { (xx) ; x e L " } - Thus 5(Xl) n 5( ^ } = { (x%) ; x e M^i-^-)}- 

There are n 2 indices j such that dim(ker(Xi — Xj)) = 1. For such an index j, 
S'(Xi) f]S(Xj) is an eigenspace of M associated to some Aj € ct(M) and conse- 
quently, Aj e i. D 

Remark 1. i) When the non-generic Eq {3J) /ias a finite number v of solutions, 

then v < ( ")• T/ie number of solutions ( n ) is known from Sylvester (1885, cf. 

PH p. 272-277];. 

m) J/ M is non- derogatory, then M admits a finite number of invariant subspaces 

and the associated Eq (0) has a finite number (eventually 0) of solutions (cf. [H 

ch. 17.8];. 

Hi) In [10] . the authors study the equation 

(7) AX 2 + BX + C = n 

and, using a method adapted from the analysis of the Riccati equation, they construct 
solvents of Eq H\). For our purpose, Eq FN is equivalent to X 2 + B' X + C = 0, 
that is Eq (0). The conclusion is very different if we consider the following non- 
unilateral equation, whose LHS linearly depends on the indeterminates, 

(8) X 2 A + BX + C = n . 

Indeed, for us, Eq $EJ) is equivalent to X 2 + BXB 1 + C" = 0, that is Eq Q), equation 
whose LHS is not a linear function of the indeterminates (cf. Section 4)- 

Definition. We consider an Eq {J) that has a finite number v of solutions. Assume 
that finding the entries of each solution of Eq fJJj is equivalent to find the roots of 
a polynomial whose Galois group over K is G. Then we denote by the solvability 
complexity (SC) of Eq l]\), the couple (u,G). 

Corollary 1. The SC of Eq 0j is (( 2 ^),S 2n )- 

Proof. It is a consequence of Theorem [T]i) and of the proof of Theorem \l\ii). □ 
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3. TWO SPECIALIZATIONS 

• Putting A — I n , B\ = B-2 = B, we specialize Eq © 

(9) X 2 + BX + XB + C = 0„. 

where the indeterminates are B, C . 

Proposition 1. Eq |3)) admits 2 n distinct solutions in the unknown X G M. n (K) 
and if n < 4, then Eq |P|) is solvable. 

Proof. Eq © is equivalent to (X + B) 2 = B 2 - C. Since C is generic, B 2 - C 
is generic and diagonalizable with non-zero distinct eigenvalues (Aj)j. There is an 
invertible matrix P such that B 2 — C = P diag(Ai, • • • , Azi)P _1 . Since X + B and 
B 2 — C commute, X = —B + P diag(jUi, • • • , /i4)P _1 where ji 2 = \. If n < 4, then 
we can explicitly calculate the (A;); and Eq © is solvable. □ 

• The equation 

(10) (XU + V)(RX + S) = n , 

where the indeterminates are the nx n matrices U, V, R, S, is a specialization of Eq 
@ with the condition CB2 1 A = B\. This equation is studied in [IB] . 
Put Y = RX + S. Since R, U are generic, they are invertible and X = R^ 1 (Y — S). 
We deduce easily that (UY) 2 + (-US + URVU~ l ){UY) = 0. The problem is 
reduced to solving the equation in Z e M n (K) 

(11) Z 2 +TZ = n 

where T is an indeterminate. Note that Eq (lll[) plays an important role when 

solving Eq © (cf. 14 ]). Clearly there are exactly 2™ solutions that commute with 

T. Yet, we go to see that there are many other solutions. 

Notation i) Let d be the Hilbert dimension of the ideal generated by the solutions 

of Eq (fTu| or of the associated Eq (fTT| ; let S be the component of the set of 

solutions of Eq (jlOllllI) that has dimension d. 

ii) Let G(r, n) be the set of r-dimensional subspaces of K n (Grassmannian) ; it is 

a homogeneous space of dimension 5 r — r(n — 7'). 

Proposition 2. i) If n is even (n = 2p), then d — S p and S consists of ( p ) affine 

varieties, whose equations have coefficients in Q[(A;)j]. 

If n is odd (n = 2p + 1), then d = S p and S consists of 2( v J affine varieties, 

whose equations have coefficients in Q[(Aj)j]. 

ii) If n < 4, then Eq ill]) is solvable. 

Proof, i) It is sufficient to consider Eq (fTTj) . Since T is generic, we may assume 

T = diag((Ai)i) where the (A-j), are indeterminates. Here M = I „" " 

Let (ej)i<2n be the canonical basis. Since T is generic, M is diagonalizable and 
ker(M) = [ei, ••• ,e n ] is the eigenspace associated to the zero eigenvalue. An 

eigenvector associated to the eigenvalue A.; is Vi = -r—ei + e n+ i. To construct a 

Aj 
Af-invariant n-vector space, we choose r < n and n — r vectors Vi x , • • • , Vi n _ r that 

generate a subspace P„_ r . Any subspace G r of [ei, • • • , e„] of dimension r does 

not intersect P n _ r and P n _ r © G r is convenient. For the sake of simplicity, we 

take the vectors V\, ■ ■ ■ , V n - r . According to Section 2, the associated solutions are 

VU^ 1 = I n n ) wnere D — diag(— Ai, • • • , — A„_ r ) and Y is an arbitrary 

\fjr,n — r ^r,rj 

(n — r) x r matrix. Each choice of G r gives birth to a solution and, when P n _ r is 
fixed, the set of solutions consists of an affine variety of dimension S r . Now sup r (<5 r ) 
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is S p if n is even and is S p = 5 p +\ if n is odd. Since there are ( n ) choices for F n _ r , 

we deduce the required results. 

ii) If n < 4, then we can explicitly calculate the (Aj)j. D 

Remark 2. The previous equation is in the form, Z 2 + BZ + C = n where BC = 
CB and B 2 — AC admits exactly 4 square roots. According to [101 p. 499], the 
solutions would be given by the usual formula for the roots of a scalar quadratic, 
that is, there would be exactly 4 solutions. In fact v is infinite and there are many 
other solutions that do not commute with T. 



4. A MORE GENERAL EQUATION 

In this section, we show a general result that we shall need in the sequel. We 
consider Eq ([T|) ; it is equivalent to a system of n 2 algebraic equations in the n 2 
unknowns (x k .i) k j 

for every i,j <n , fij(X) = 0. 

Definition. Let Xq be an isolated solution of Eq flj). Let Ox be the local ring in 
Xq and {{fi,j)i,j)x be the ideal generated by the {fi,j)i,j in Ox - The multiplicity 
of X is 

the dimension of Ox /{{fi,j)i,j)x , as a K-vector space. 

Example. Consider the solution (0,0) of the system in C 2 : {x 2 = 0,y 3 = 0} ; 
the quotient C[x, J/](o,o)/(* >ir)(o,o) admits the basis {l,x,y,y 2 ,xy,xy 2 } and has 
dimension 6. For the system {y 2 — x 5 — 0,x 2 — y 5 = 0}, the multiplicity of the 
solution (0, 0) is 4 because it is locally equivalent to {x 2 — 0, y 2 = 0} . 
If we know a Grobner basis of the ideal generated by the {fi,j)i,j, then fortunately, 
we can calculate (Maple does it) the total number of solutions of Eq {IJ) with mul- 
tiplicity. For our second system, this number is 25, that is 4 solutions in (0, 0) and 
the 21 solutions of the subsystem {y 21 — l,x = y 13 }. 

Theorem 2. The n x n matrices (A%, ■ ■ ■ ,A k ,B) are indeterminates and g is a 
non-zero polynomial over K in fe + 1 non-commuting indeterminates such that each 
monomial of g contains explicitly the variable X. We consider the equation 

(12) g{M,- ■ ■ ,A k ,X) + B = 0„ in the unknown X e M n {K). 

Then Eq 112(1 admits no singular solutions, that is, there are a finite number (even- 
tually 0) of solutions, and all have multiplicity 1. 

Proof. Assume that there are an infinity of solutions or there is an isolated multiple 
solution ; according to the intersection theory (cf. |5],[7]), that is equivalent to the 
existence of a singular solution Xq. Thus Xq satisfies 

det(||(A 1 ,...,A fc ,Xo))=0, g(A 1) >>-,A k ,X ) = -B. 

The hypersurfaces (gij(Ai,- ■ ■ ,A k ,Xo) = —bij)ij have no common components 
and, according to Bezout's Theorem (cf. [5]), the equation g(A\, ■■ ■ ,Ak,X ) = —B 
has a finite number of solutions in Xq that depend on B. Thus, if such a solution 

dg 
Xq satisfies det(jrr^ (Ai , • • • ,A k , Xoj) = 0, then we deduce a relation linking the 

oX 
(Ai)i and B, that is a contradiction. □ 
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5. Random specializations and formal calculus 

In this section and also in Section 7, we show the interest of the formal calculus 
in the study of the generic matrix equations. Moreover the notion of "thin set" is 
usefull to construct, in a probabilistic way, a generic matrix equation. 
• Tests on Eq @. Let n € [2,4]. Of course we can calculate explicitly xm- 
Practically, since we cannot explicitly calculate cr(M) (except when n = 2), what 
can we do after ? A possible answer is using the Grobner basis theory (cf. [3]). 
According to Section 2, it is sufficient to consider Eq ((3j) . Since B is generic, we may 
assume that B — diag((Aj)j), where the (Aj)j are distinct and can be calculated by 
radicals. Let C — [ctj]. Thus we must solve a system of n 2 algebraic equations 
in the n 2 unknowns (a?t,j)t,j and in the the n 2 + n indeterminates (A,),-, {ci,j)i,j- 
Let L be the field Q((Aj),, {ci,j)i,j) and t be the trace of X. Since there are a 
finite number of solutions in X, then the Hilbert dimension of the ideal generated 
by these polynomials is equal to n 2 + n, the number of indeterminates. We use 
the package "Groebner" of Maple. When the calculation is complete, we obtain a 
Grobner basis that triangularizes the system as follows 

P(t) = ; for every i, j , x lJ U lJ (t) = V i; j(t), 

where P, a polynomial over L, has degree ( 2 ™) and S2n as Galois group and the 
(Uij), (Vi.j) are polynomials over L of degree ( 2 ™) — 1 such that, for every i,j, Uij 
and P have no common roots. Unfortunately the calculation is complete only when 
n = 2. When n = 3, too much memory is required, even with a cluster. Thus, 
when n = 3, 4 we must specialize the matrices B,C. 

The reason that we choose the unknown t is the following : the minimal polyno- 
mial over L of t is in the form P(z) = z T +Xh=o a i z% where t = ( ") . At least when 
n < 4, Q(u) — P(u — ) is an even polynomial. Therefore Q(u) = R(u 2 ) where 

T 

R is a polynomial of degree r/2 and the problem is reduced to solving a polynomial 
of degree 4( ")• In particular, when n = 2, R has degree 3 and we obtain easily 
the values of t. When n = 3, we consider this random instance of specialization 



B= I I , C 
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Here R is the following irreducible polynomial of degree | (T) = 10 
R(z) = z w - I0z 9 + 57z 8 - 263z 7 - 2350z 6 - 5817z 5 

+31627z 4 + 289793z 3 + 373833z 2 - 600264z + 60516. 

Let H be the Galois group of R over Q. Using the software Magma and according 
to the notation of [3], we obtain H w "10 T 32". This group is isomorphic to 56- 
The Galois group of P is also S&, as for the generic system, that is not surprising 
for the reasons below. Note that, if n — 4, then the calculation is complete only if 
we randomly choose for B a diagonal matrix. 

Recall that the specialization of the indeterminates results in a Galois group 
which is a subgroup of the original Galois group. Yet, if the choice is random, then 
(except if we are unlucky) we obtain the total group. To study this fact, we need 
the concept of "thin set", in the sense of J. P. Serre (cf. [TTJ p. 19]) ; its precise 
definition is useless for our purpose. Write £ € -Pi(Q) = QU{oo} as the equivalence 
class of (x, y) € 1? \ {(0, 0)} with gcd(x, y) = 1. This can be done in a unique way, 
up to a choice of sign. The height of £ is defined to be height (£) = sup(|x|, \y\). 
The following shows that when we randomly choose a point in P\ (Q) , it is almost 
never in a fixed thin set 
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Proposition. 17, p. 26] i) The number of points in Pj.(Q) with height less than 

N is asymptotically -p-A 2 . 

it) If A a Pi(Q) is a thin set, then the number of points of A with height less than 

N is inO(VN). 

We consider a polynomial, the coefficients of which depend on indeterminates. 
The following result shows that a random specialization of the indeterminates al- 
lows to retrieve the total Galois group of the polynomial. Let V be an algebraic 
irreducible variety over Q. Let V(Q) denotes the Q-rational points of V and Q(V) 
denotes the function field of V. 

Theorem. Q21 p. 24] Let f(x) = x n + a n _ia;" -1 H \-a a , where a, e Q(F), be 

an irreducible polynomial over Q(V) that has G as Galois group. Then there exists 

a thin set A C V(Q) such that if t ^ A, then 

i) t is not a pole of any of the (a{). 

ii) ft(x) = x n + a n -\{t)x n ~ l + ■ • • + ao(i) is irreducible over Q. 

Hi) The Galois group of ft is G. 

If the parameter N is not a large number, then the probability of falling into a 
fixed thin set is not negligible, as we see about the following test where the entries 
are randomly chosen in [— 2, 2J. 

Remark 3. When n — S, there exist solvable instances {B,C) of Eq ^Bj), such that 
BC — CB is not nilpotent (cf. Theorem^ when B,C commute), as this one 

/ 2 -12^ 
B = diag(l,-1,2) , C = -1 1 

\1 2 y 

Here, R = R^Rq, where i?4,i?6 o^e irreducible polynomials of degrees 4,6 and Rq 
is solvable. 

• Tests on Eq (gj. Now we consider Eq Q : X 2 + BxXB 2 + C = 0„, in the 
unknown X E A4 n (K), where the n x n matrices B\,Bi,C are indeterminates. 
Unlike to Eq ([2]), its LHS does not linearly depend on the indeterminates and, 
consequently, we do not know any explicit form of the solutions of Eq Q ; thus we 
study this equation using formal calculus for n — 2,3. We obtain that Eq (QJ has 
a larger SC than that of Eq ©. 

Proposition 3. For any n > 2, Eq (CI) has a finite number of solutions and is 
non-solvable. 

Proof. According to Theorem [21 Eq (JU) has a finite number of solutions ; thus, 
according to the Grobner basis theory, there is (i, j) and a polynomial p, with 
coefficients in K, such that the solutions Xij are the roots of p. If a specialization 
of Eq (UJ) admits a finite number of solutions Xq, then the solutions Xo-ij are the 
roots of po, the specialization of p. Since the number of solutions X is finite, p is 
not identically zero. Suppose that Eq Q is solvable ; therefore the polynomial p is 
solvable and, consequently, po too. 

• n = 2. Using Maple, we consider a specialization -Bi,o,-B 2i o,Co of Eq (UJ) and 
construct a Grobner basis of the ideal generated by the associated 4 algebraic 
equations in 4 unknowns. The polynomial po, associated to x±.i, has degree 8 and 
its Galois group is Sg, that is a contradiction. 

• n > 2. Eq ([3]) is a specialization of Eq (g]). According to Theorem Q] Eq ^ has 
a finite number of solutions and is non-solvable ; thus Eq (J4J is non-solvable. □ 

More precisely we do specializations as follows : the entries of B\,B?,,C are 
randomly chosen and we calculate the polynomial po associated to X\\. Almost 
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all the tests are such that po has degree 8 (resp. 56) and Sg, (resp. S^) as Galois 
group over Q. Thus, according to the Serre theorem, we conjecture the following 

Conjecture 1. When n = 2 (resp. n — 3), the SC of Eq Q) is (8,Sg) (resp. 
(56,S 56 )j. 

6. When the indeterminates commute 

Now we assume that the indeterminats commute. Let Bq, ■ ■ ■ ,-Bfc-i be n x 
n commuting matrices that are indeterminates and K = Q(Bq,--- ,Bk-i). We 
consider the following unilateral matrix equations in the unknown X £ M. n (K) 

(13) X k + X k - 1 B k ^ x + ■■■ + XB X +B Q = 0„. 

(14) X k + B^X^ 1 + ■■■ + B X X + B a = 0„. 

Note that, in the previous equations, the LHS is a linear function of the indetermi- 
nates and the monomials of type BX k and X l C are not mixed. 

Theorem 3. We consider Eq H3\) or Eq TTM where the (Bi)i are nxn commuting 
matrices that are indeterminates. Then there are k n solutions and any solution 
X £ Ai n (K) commutes with the (B{)i. 

Proof. By transposition, solving Eq (J14I) is reduced to solving Eq (fTB")) . Thus we 
study only Eq (|1 3[) . Since Bq is a generic matrix, it is similar to a diagonal matrix 
diag(Ai, • • • , A„) whose distinct eigenvalues are conjugate over K and the commu- 
tant of Bq is if [Bo]- Thus, for every j < k—1, there exists a unique polynomial Pj 
of degree n — 1 and with coefficients in K such that Bj = Pj(Bq). Let (ej)i<„ be 
a basis consisting of eigenvectors of Bq. Thus for every i < n 

(X k + Pu^iX^- 1 + ■■■ + Pi(\i)X + Xil n )ei = 0. 

Put 

6{x, X) = x k + Pfc-itA)^- 1 + • • • + Pi (X)x + X 

and, for every i, let ^ i, • • • , /ii j. be the roots in the unknown x of 8(x, Xi). 

If the (fii,j)i< n ,j<k are not distinct, then there is i < n such that the discriminant 

discvim(9(x 1 Xi),x) = 4>{Xi) is zero, 

or there are p < q < n such that the resultant 

rcsult(6'(a;, A p ), 9{x, X q ), x) — ip{X p , X q ) is zero. 

Note that <f>, ip are polynomials with coefficients in K that depend only on the 
(Pj). We consider symmetrizations of (j)(Xi), n CTe 5 n 0(A CT (j-)) = 0, or of ijj(Xp,X q ), 
n cre 5 n , 0(A cr (p), A CT ( 9 j) = 0, that are polynomial relations linking the matrix Bq and 
the polynomials (Pj). We obtain a contradiction if we show that <fi(Xi) and ip(\i \) 
are not identically zero. 

We specialize the (Bi)i into the (B'^i £ Ai n (Q) , putting, for every r < n, 
X r = r and, for every < s < k — 1, P s (x) — x s+1 mod XB a ( x )- Then a 
calculation shows that the associated (/ij ~)i< n ,j<k are distinct. Note that every 

fe-tuple of commuting matrices in a neighborhood of (B'^i in A / l n (Q) satisfies the 
previous property. 

We return to the case where the (Bi)i are commuting indeterminates. There 
are distinct /i; : i, • • • , /i^fe such that (X — fii t il n ) ■ ■ ■ (X — / Ui l fc/„)ej = 0. Thus there 
is ji < k such that [i^ £ cr(X). Since the (Hi,j)i< n ,j<k ar e distinct, X has the 
following simple eigenvalues : {(Hj^iKn- We deduce that, for every i, Xei = Mi,j« e « 
and X commute with the (Bi). It remains to solve n polynomials of degree k. □ 
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Remark 4. An essential argument in the previous proof is that the {t L i,j)i<n,j<k 
are distinct. Thus we may replace, in Eq U3\) . the polynomial X k + X k ~ 1 Bk-i + 
• • • + XB\ + Bq with a sparse polynomial. Yet, necessarily, B$ must appear (cf. Eq 
Ml}) in Section 2) ; else x = is a common root of the polynomials (9(x, \i))i- 

Corollary 2. Consider the Riccati Eq $£) where the indeterminates A,B\,B 2 ,C 
are commuting n x n matrices. Then there are 2™ solutions, in X € Ai n (K), and 
any solution commutes with the indeterminates. 

Proof. Using Section 3, a solution X of Eq ([2]) can be write X — A~~ l Z — A~ l B 2 
where Z is a solution of the equation Z 2 + DZ+E = 0„ with D = —B 2 +AB\A~ 1 = 
-B 2 + B U E = -AB 1 A- 1 B 2 + AC = -B X B 2 + AC. The indeterminates D, E are 
only linked by the relations of commutativity : DE — ED = 0„. According to 
Theorem [3J Z is a polynomial in E. Moreover A,A ,Bi,B% are polynomials in 
C. Then Z is a polynomial in C and X too. □ 

Consider, for n = 3, the equation of degree 4 

(15) X 4 + AiXA 2 XA 3 X + X 3 + X 2 + XA A XA 5 + A 6 XA 7 + B = 3 

where the 3x3 matrices (Ai)i,B are commuting indeterminates. Using Maple, we 
specialize Eq (|T5|) and we obtain almost always that the solutions inle M.^(K) 
commute with (Ai)i,B. Other numerical experiments do seem to indicate that the 
following surprising result is true 

Conjecture 2. Let g be a non-zero polynomial over K in k + 1 non-commuting 

indeterminates (u\, • • • , Ufc+i) such that 

i) g(ui,-- ■ ,Uk,0) =0 

and for every i < k 

ii) Ui appears exactly one time in the expression g{u\, ■ ■ ■ ,Uk+i)- 

Hi) The degree of g with respect to Ui is 1. 

We consider the equation 

(16) g{Ai, ■ ■ ■ ,A k ,X) + B = O n in the unknown X e M n (K) 

where the n x n matrices A\,--- ,A^,B are commuting indeterminates, that is, 
we assume that the only ones relations that link the indeterminates are those of 
commutativity. Then any solution of Eq il6}) commutes with the indeterminates. 

According to Conjecture [2j any solution of the commuting Eq (j4]) would com- 
mute with B\,B 2 ,C ; in fact, we can prove this result, using Maple, for n < 3. Yet, 
the following shows that Conjecture [5] does not work if condition ii) or Hi) is not 
fulfilled. 

Proposition 4. We consider the equation : X 2 + BXB + C = 2 (that is a special- 
ization of Eq (Opj where B,C are 2x2 commuting matrices that are indeterminates. 
Let L = K(a(B), o~(C)). There are four trivial isolated solutions that commute with 
B,C. Moreover there are infinity many solutions X £ Ai 2 (L(u)) that depend on a 
parameter u ; these solutions do not commute with B, C. 

Proof. We may assume that B, C are diagonal matrices. Using Maple, we construct 
a Grobner basis of the ideal generated by the associated 4 algebraic equations in 4 
unknowns and 4 indeterminates. Since the Hilbert dimension of the ideal is 5, the 
non-trivial solutions depend on 5 — 4 = 1 parameter. □ 

Remark 5. More generally, the previous equation admits, in A4 n (K), non-trivial 
solutions that depend on n — 1 parameters. 
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7. Quadratic equation in matrices 

The general quadratic matrix equation, in the unknown X = [Xi t j] € M. n {K), 
has the form (cf. [15] ) 

r 

(17) ^(AiXBiXd + DiXEi) + F = 0„ 

i=l 

where the n x n matrices (^4»), (-Bt)j (Cj), (Dj), (-Ej), i* 1 are indeterminates. 

Assume n = 2. Eq (fT7]l is equivalent to 4 algebraic equations of degree 2 in 
the 4 unknowns (xi.j)ij. Using a supplementary unknown T, we homogenize the 
previous equations and we seek the intersections of 4 quadric hypersurfaces (Hi)i 
in the projective space P±{K). Let v be the number of solutions of Eq ([T7|l in 
P4(i4T), counted with multiplicity, and by considering points at infinity. According 
to Bezout's Theorem, 

i) or v is finite and v = 2 4 = 16. Moreover, if we know 16 isolated solutions, then, 
necessarily v is finite and there are no other solutions. 

ii) either v is infinite and the hypersurfaces (Hi)i have a common algebraic compo- 
nent. Then the sum of the degrees of the irreducible components of their intersection 
is at most 16. 

Note that the solutions such that T = satisfy the homogeneous equation 

r 

(18) ^(AtXBiXd) = 2 . 

1=1 

Thus Eq (|17p admits no solutions at infinity when Eq (|18l) admits the unique 
solution X = 2 . 

Example 1. The set of solutions of the Riccati Eq |Q|) consists of 6 isolated solu- 
tions in K 4 and, at infinity, the set {X ^ | XAX = 2 } = {X ^ | (XA) 2 = 2 } 
; the previous set is a blunted cone of dimension 2, that is a curve in the hyperplane 
T = QofP 4 (K). 

Example. Consider the quadratic equation AXB1X+XB2X+DX + F = 2 where 
the 2x2 matrices A, B\,Bi,D,F are indeterminates. Numerical experiments seem 
to indicate that this equation admits 8 distinct solutions in K . At infinity, we study 
the equation ip(X) = AXB\X + XB-^X = 2 . We obtain almost always 4 distinct 

solutions in P^(K) ; if Xq is such a solution, then -j— (Xq)(Xq) = 2^(Ao) = 2 . 

oX 

dtp 
Thus rank(— — (Xq)) < 4 and Xq is a multiple solution ; in order to obtain 16 

oX 

solutions in P±(K), its multiplicity would be 2. 

We give an instance of a generic quadratic matrix equation that admits 16 iso- 
lated solutions in K . 
Let Kd denotes the set of polynomials of degree d with cofficients in K. 

Theorem 4. The 2x2 matrices A, B\, B2,C, D, F are indeterminates. We con- 
sider the equation in the unknown X G M.i{K) 

(19) AXB 1 X + XB 2 X + X 2 C + DX + F = 2 . 

The SC of Eq (Wj) is (16,S*i 6 ). 

Proof. • We specialize Eq (fi"9|) and use the "Groebner" package of Maple. If we 
randomly choose the indeterminates, then we obtain almost always the following 
i) A XBi_ X + XB 2: qX + X 2 Co = 2 implies X — 2 , that is our specialization 
has no solutions at infinity. 
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ii) a Grobner basis associated to Eq (|19p. in the unknown Xq = [#o;j,.j]j has the 
form 

Po(»0;i,i) = 0, for every (r, s) 7^ (1,1), £o ; r, s = Po;r,*(#0;i,i) 
where po € K\& has Si 6 as Galois group and po-, r ,s G if 15- 

Thus we obtain 16 distinct solutions in K . 

• According to Theorem^ Eq (|19p has a finite number of solutions in M.z{K ). As 
in the proof of Proposition [21 there is (i, j) and a polynomial p, with coefficients in 
if and deg(p) < 16, such that the solutions Xij are the roots of p. If a specialization 
of Eq Q admits a finite number of solutions Xo, then the solutions Xo-ij are the 
roots of Pq. Since the number of solutions Xq is finite, po is not identically zero 
and deg(po) < deg(p). Our experiments show that deg(p) = 16 and that the Galois 
group of p is Si6- Finally Eq ([15]) has 16 distinct solutions in M 2 (K) and each 
value of Xij gives birth to a unique solution X . □ 

The nx n matrices A, Bi,B 2 ,C are indeterminates. We consider the homoge- 
neous equation, in the unknown X e •Mn(.K'), which is associated to Eq (|19D 

(20) AXSiX + XS 2 X + X 2 C = 0„. 



Corollary 3. If n = 2, then Eq H20\) has only the zero solution. 

Proof. According to Theorem UJ Eq (|19p has no solutions at infinity. Therefore, 
the homogeneous Eq ([20)1 has no non-zero solutions. □ 



Assume n = 3. We look for a generic quadratic equation that is maximal in the 
sense that it has 2 9 = 512 isolated solutions in M.^(K). Unfortunately Eq (|T9|) 

is not convenient. The set of non-zero solutions of Eq ([2"u]) is a cone in K \ {0}. 
Therefore we consider solutions in Pg (K) . When we specialize Eq (|20l) , the Hilbert 
dimension of the associated ideal is almost always 1, that is, there is at least one 
curve of solutions in Ps(K ). To obtain a maximal equation, we proceed as follows 
: first, we seek a homogeneous equation ip(X) = O3 that has no non-zero solutions 
; then, we consider the equation ip(X) + DX + F = O3 where the 3x3 matrices 
D,F are indeterminates. According to Theorem^ the previous equation has only 
isolated simple solutions ; thus it has 512 distinct solutions in A4s(K). According 
to numerical experiments, we conjecture the following 

Conjecture 3. The 3x3 matrices (Ai), (Bj), (Ck), D, F are indeterminates. Then 
the equation in the unknown X £ M.?,{K) 

ip(X) = AyXBxX + A 2 XB 2 X + XB 3 Xd + XB A XC 2 + XB b X = 3 

has no non-zero solutions, or equivalently, the solution X = O3 has multiplicity 512. 
Moreover the SC of the equation ij)(X) + DX + F = O3 is (512, S512). 

8. Instances of Riccati equations 

In this section, we adopt another point of view. Let n — 2 and v denote the 
number of solutions, with multiplicity and not at infinity, of Eq ©. Here the 
matrices A,Bi,B 2 ,C are numeric and fixed in M 2 (C) ; the only condition that is 
required on these matrices is A ^ Q 2 . We look at the different forms of the set of 
solutions of a non-generic Riccati equation. The results that follow are known by 
V.V. Palin but it seems that they were not published. 

• We know that if v is finite, then v < 6. For every r < 6, we give an instance of 
Eq @ such that v = r. Recall that if M is non-derogatory, then v is finite. 

r = 0. X 2 — I I = O2 or XAX + C = 2 where A is not invertible and C is 

invertible. 
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r = l. x(l l\x+(~Q ^ X +(~i -l) = ° 2 - ^W = {-1.-1. 0,0} and 

M is non-derogatory. The unique convenient M -invariant subspace of dimension 2 

is ker(A/ 2 ). 

Question : Does such a simple solution exist when A is invertible ? 

, 2 . „ . (\ 



2. X 2 + X + I 1 = 2 . cr(A/) = {A, A, A, A} and M is non-derogatory. We 



obtain 2 solutions : Xo,Xq. 

r = 3. X 2 + I J A+ I J = O2. The unique solution is X = I ' 

and has multiplicity 3. a(M) = {1, 1, 1, —2} and M is non-derogatory. If we perturb 

the entries of the coefficients of X and X°, by adding small complex numbers, 

then according to Section 2, we obtain (in general) a generic equation that has 6 

isolated solutions ; 3 solutions are in a neighborhood of Xq and 3 solutions are in 

a neighborhood of infinity, that is they have at least one large entry. 

r = 4. Use a generic commuting Eq (J2]) (cf. Theorem [3]). M has four distinct 

eigenvalues. 

r = 5. X 2 +( n 1 I X + I _ 1 I = 02- M has four distinct eigenvalues. 

r = 6. Use a generic Eq (J2J) (cf. Theorem [T). M has four distinct eigenvalues. 
r — 6. (communicated to us by V.V. Palin). 

*(! !)*-(§ i)^(§i)-^ 

M is nilpotent and non-derogatory ker(M 2 ) is the unique 2-dimcnsional M- 
invariant subspace. The unique solution is X = O2 and has multiplicity 6. 

Remark 6. More generally, consider the equation XAX + BX — XB = O2 where 
the 2x2 matrices A, B are indeterminates. Then we can show that the previous 
equation admits the trivial solution X = O2 with multiplicity 4 and 2 supplementary 
non-zero solutions. 

• Assume that v is infinite. As we shall see, the Hilbert dimension d of the ideal 
generated by the solutions is 1 or 2. 

Remark 7. In Section 3, we saw that, if A is invertible, then Eq @) can be 
rewritten in the form of Eq $) : X 2 + BX + C = 2 . One has dct(M - XI 2 ) = 
det(A 2 /2 — AS + C). Assume that B, C commute ; then we may assume that B,C 
are upper-triangular and there are orderings (fix, ^2) and (y\, v-i) of o~(B) and o~(C) 
such that the eigenvalues of M are the roots of the polynomial (A 2 — fi\\ + ^i)(A 2 — 
fi2^ + V-i)- Since v is infinite, the previous polynomial has a double root. 

d = 1. Use a generic Eq (II) (cf. Proposition [2). a(M) = a(T) U {0,0} (3 
distinct eigenvalues). M is derogatory and diagonalizable. S consists of 2 non- 
intersecting straight lines in the direction of 2 nilpotent matrices. 
Question : Does there exist an instance of Eq @ such that d = 1 and that the 
variety of solutions is not flat ? 

d = 2. X 2 - I 2 = 2 . The set of solutions is {I 2 , -J 2 } U{ ( ° b ) \ a 2 + be = 1}, 

that is 2 isolated points and a hyperboloid of one sheet of dimension 2. a(M) = 
{A, A, /1, fi} where A^/i. M is diagonalizable and derogatory. 

d = 2. XAX = O2 where A 7^ O2. M is nilpotent and derogatory. We obtain a 
cone of dimension 2 when A is invertible (see Example [I) and else, considering the 

cases A = I I and A = n ), we obtain two intersecting planes. 
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The previous instance shows that the set of solutions at infinity is always a curve 
in P4, that implies that the maximal Hilbert dimension of the ideal generated by 
the solutions, not at infinity, is 2. We deduce the following 

Proposition 5. We consider Riccati Eq (0) where n = 2 and A € A^C) \ 

{O2}, Bi, B2, C € A^2(C) are given fixed numeric matrices. There exist instances of 

the previous matrices such that the set of solutions, not at infinity, of the associated 

Eq (0) is one of the following 

i) 0, 1, • • • or 6 elements. 

ii) A variety of Hilbert dimension 1. 

Hi) A variety of Hilbert dimension 2. 

Conversely, each instance of Eq (0) has a set of solutions that has one of the 

previous forms. 
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